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ABSTRACT 

;  v 

V  The  ratio  of  the  "a"  phase  and  "max"  phase  of  presumed  Shagan  River  contained 
and  cratering  explosions  are  studied  across  the  WWSSN  network  of  short  period  sta¬ 
tions.  The  amplitude  of  the '"a"  phase  of  the  presumed  cratering  explosion  of  January*2- 
15,  1965  is  found  to  be  systematically  larger  in  comparison  to  the  "max"  amplitude  of 
this  cratering  explosion  when  compared  to  contained  explosions  in  the  vicinity  of  the 
cratering  explosion^Tn--other  words^the  log  of  the  ratio  of  the  amplitudes  of  thc'-"a"-'knd 
"max"  phase,  log(Pmax/Pa),  is  on  average  smaller  for  the  Janbary^  15,  1965  cratering 
explosion  than  for  any  of  the  contained  explosions  studied.  Several  methods  were  used 
to  study  the  systematics  of  the  log(Pmax/Pa)  value  across  the  WWSSN  network.  The 


r  t  W  '  r  , 

preferred  method  for  ^determination— of^the  average  log(Pmax/Pa)  was  a  maximum- 
likelihood  method  that  includes  the  effects  of  data  truncation  of  clipping  for  Pmax  and 
non-detection  of  Pa.  Assuming  that  the  phase  amplitude  of  the  cratering  explosion 
is  unaffected  by  the  influence  of  the  non-linear  free  surface  interaction,  the  corrected 
magnitude  for  the  cratering  explosion  is  0.17  to  0.27  magnitude  units  larger  than  its  5.88 
WWSSN  network  magnitude.  If  the  cratering  explosion  is  credited  with  a  yield  of  125 
KT  then  a  125  KT  contained  explosion  in  the  Shagan  River  test  site  should  produce  a 

b  _ _ 

WWSSN  network  (m%\  between  6.05  to  6.15.  Assuming  a  log(yield)  slope  of  1.0  the 
analysis  predicts  a  150  KT  explosion  should  produce  a  WWSSN  network  m»  between 
6.13  and  6.23.  The  analysis  assumes  that  the  Pa  phase  amplitude  is  the  same  for  crater¬ 
ing  and  contained  explosions  of  the  same  yield,  and  the  WWSSN  network  at  distances 
between  20  and  90  degrees. 
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INTRODUCTION 

A  detailed  determination  of  cratering  and  non-cratering  short  period  P-wave  magni¬ 
tudes  was  conducted,  in  order  to  better  understand  the  short  period  magnitude  biases 
between  different  test  sites.  The  P-wave  mb’s  were  determined  using  the  initial  "a" 
phase,  the  "b"  phase  and  the  maximum  P-wave  amplitudes  in  the  first  5  seconds.  Mag¬ 
nitudes  based  on  amplitude,  A,  and  amplitude  corrected  for  period,  A/T,  were  con¬ 
sidered.  The  differential  between  these  different  methods  of  P-wave  amplitude  determi¬ 
nation  was  studied  as  a  network  average  and  on  a  station-by-station  average.  It  has 
been  suggested  that  the  initial  displacement  of  the  tclescismic  P  wave  is  not  as  affected 
by  the  free  surface  effects  that  are  manifested  in  the  pP  arrival.  The  detonation  of  a 
shallow  explosion  results  in  a  non-linear  "reflection"  of  energy  at  the  free  surface  that 
must  appear  near  the  time  of  the  linear-elastic  pP  wave.  The  scaling  of  this  pP  with 
depth  and  yield  remains  a  major  problem  in  the  understanding  of  the  tclescismic  short- 
period  waveforms.  The  P-fpP  waveform  is  certainly  observed  for  deep  explosions  while 
shallow  explosions  show  less  definitive  evidence  of  the  P+pP  waveform.  It  is  argued 
that  the  "a"  phase  is  less  affected  by  the  free  surface  reflection  even  for  cratering  explo¬ 
sions  where  there  is  no  true  reflection  of  energy  at  the  free  surface.  This  argument 
implies  that  the  "a"  phase  will  scale  more  linearly  with  yield  than  the  "b"  and  maximum 
phases  commonly  used  in  magnitude  determination.  The  depth  correction  for  the  "a" 
phase  is  therefore  unnecessary,  while  it  may  be  considered  necessary  for  "b"  and  "max" 
phases. 

Given  that  the  Pmax/Pa  amplitude  ratio  can  be  determined  for  the  "normally" 
buried  and  cratering  explosion  created  teleseismic  P  waves,  a  check  on  the  yield  of  the 
East  Kazakh  explosions  can  be  formulated  from  the  difference  between  the 
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!og(Pmax/Pa)  for  the  scale-depth  explosion  and  the  cratering  (125  KT)  explosion  at 
East  Kazakh.  A  bias  for  the  cratering  explosion  as  measured  using  the  log(Pmax/T)  or 
log(Pmax)  as  a  measure  of  magnitude  can  be  determined.  We  can  define  this  predicted 
"bias"  as  5,  where 

S  -  log  (Pmaz  /  P.i  )„.<,  _*,«*  -  log  (Pmaz  /  Pa 
The  mb  magnitude  of  a  scale-depth  explosion  of  the  same  yield  as  the  cratering  explo¬ 
sion  of  Jan.  15,  1965  would  then  be 

m^iUKT  “  m^trmt ir  +  5. 

Teleaeismic  P  waves  were  synthesized  for  the  Von-Seggern  Blandford  (1972)  source 
time  function  observed  through  a  WWSSN  short  period  instrument.  The  pP-P  delay 
time,  r  in  seconds,  was  assumed  to  scale  by  either  r  t=0.07  W'/  *  or  r  j=0.12  W'1' 1  (  r  in 
seconds,  W  in  KT).  Attenuation  parameterized  by  t  *  values  of  0.20  and  0.45  seconds 
was  introduced  into  the  synthetic  waveforms.  The  pP  was  assumed  to  be  a  "linear" 
reflection  with  decreasing  amplitude  with  increasing  frequency  as  in  Blandford  et  al. 
(1984).  The  synthetics  were  also  computed  for  no  pP  contribution.  These  "linear" 
models  predict  log(Pc/Pa)  values  greater  than  the  values  for  no  pP  as  listed  below. 


- ~~i - 

t 

r 

log  {Pci  Pa)„rt_*,tt-/oy  {Pc  /  Pa)M.fP 

(sec) 

(sec) 

0.45 

0.37 

0.22 

0.45 

0.64 

0.19 

0.20 

0.37 

0.24 

0.20 

0.64 

0.14 

Depending  upon  the  choice  of  t '  and  the  scale-depth,  the  estimated  value  of  E 
varies  between  0.14  and  0.24.  The  same  calculations  were  conducted  for  the  LRSM 
short  period  instrument  which  has  a  broader  bandwidth  than  the  WWSSN  short  period 
instrument  with  virtually  the  same  results.  We  would  expect  that  pP-P  times  may  vary 
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over  a  test  site  if  the  testing  practices  or  geology  are  not  uniform.  We  would  expect  that 


the  effective  ( *  would  vary  over  a  network  of  stations. 


In  order  to  predict  the  5  for  150  KT  events  at  200  meters  (cratering)  and  500 
meters  (scale  depth),  while  beginning  to  take  into  account  the  non-linear  effects,  we 
examine  the  suite  of  granite  calculations  made  by  Dache  et  al.  (1980).  For  a  set  of 
depths  ranging  from  150  meters  to  1000  meters,  they  predict  tclcscismic  waveforms  for  a 
150  KT  event  observed  with  a  KS-3C000  seismometer  and  al’  =  0.8  sec.  The  following 
predicted  differences  between  mb(Pc)  and  mb(Pa)  for  these  calculated  depths  are  taken 
from  Figure  15  and  Table  5  of  Bache  et  al.  (1980). 


h 

mb(Pc)-mb(Pa) 

(meters) 

159.4 

0.49 

207.2 

0.50 

253.0 

0.50 

398.5 

0.57 

531.3 

0.62 

797.0 

0.69 

1000.0 

0.63 

The  difference  between  the  shallow  source,  about  200  meters,  and  the  deep  source, 
about  500  meters,  is  5  =  0.12  magnitude  units.  Consequently,  given  the  assumptions  in 
the  Bache  et  al.  (1980)  model  we  would  predict  that  the  "a"  phase  of  the  cratering 
explosion  is  relatively  larger,  compared  to  the  "c"  phase  by  about  0.12  magnitude  units 
for  a  similarly  sized  explosion  near  scale-depth.  This  is  slightly  less  than  predicted  by 
the  "linear"  model  presented  above. 

Several  factors  may  conspire  to  obscure  the  predicted  !og(Pmax/Pa)  distribution 
for  a  real  network  recording  real  explosions.  These  factors  include  the  variance  of  ampli¬ 
tude  data  due  to  focussing-defocussing,  the  variance  of  attenuation  due  to  source-receiver 
paths,  and  the  generation  of  P  coda  due  to  scattering.  P  coda  may  be  such  that  the 
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"Pmax"  in  the  first  5  seconds  of  the  P  wave  is  not  the  "c"  phase  as  in  the  case  of  syn¬ 
thetic  P  waveforms  in  a  horizontally  layered  earth.  The  following  section  attempts  to 
determine  if  the  log(Pmax/Pa)  is  "stable"  for  East  Kazakh  events  as  a  function  of  mag¬ 
nitude  and  locations  in  the  proximity  of  the  cratering  event  of  Jan.  15,  1965.  Following 
the  determination  of  the  stability  of  log(Pmax/Pa),  the  value  of  b  is  estimated  for 
events  with  mb  «  5.9  and  epicenters  near  the  January  15,  1965  cratering  event. 
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DATA  ANALYSIS  OF  TELESEISMIC  P  WAVES 

The  Jan  15,  1965  cratering  explosion  in  East  Kazakh  had  an  announced  yield  of 
125  kilotons  (Nordyke,  1973).  The  explosion  crater  (and  the  lake  it  became)  is  clearly 
visible  on  LANDSAT  images.  The  event  has  been  used  by  several  researchers  to  re¬ 
locate  East  Kazakh  events  and  calibrate  telescismic  networks  for  the  computation  of 
travel  times  from  the  East  Kazakh  test  site  (Rodean,  1979;  Shore,  1982).  Several  events 
of  similar  yield  have  occurred  within  several  kilometers  of  the  crater.  We  have  used  the 
WWSSN  network  to  determine  the  average  station  differentials  between  mb(b), 
mb(max),  and  mb(a)  for  a  number  of  East  Kazakh  events.  Several  of  these  events  are  in 
close  proximity  to  the  location  and  yield  of  the  Jan  15,  1965  cratering  explosion.  The 
ISC  origins  of  the  E.  Kazakh  events  studied  in  this  report  are  listed  in  Table  I. 


Table  I.  ISC  ORIGINS 

Event 

Origin  Time 

N.  Lat. 

E.  Long. 

mb 

S  (km) 

15jan65 

05:59:58.4 

49.88 

78.96 

5.8 

0.0 

19jun68 

05:05:57.4 

49.96 

79.05 

5.4 

11.0 

1 lsep69 

04:01:57.5 

49.77 

78.03 

5.0 

68.1 

30nov69 

03:32:57.3 

49.94 

78.98 

6.0 

6.8 

14dec73 

07:46:57.1 

50.03 

79.02 

5.8 

17.2 

27apr75 

05:36:57.2 

49.94 

79.02 

5.6 

8.0 

04jul76 

02:56:57.5 

49.85 

78.97 

5.8 

3.4 

07dec76 

04:56:57.5 

49.87 

78.89 

5.9 

5.2 

lljun78 

02:56:57.7 

49.88 

78.81 

5.7 

10.8 

15sep78 

02:36:57.3 

49.91 

78.94 

6.0 

3.6 

The  S  is  the  distance,  km,  from  each  ISC  epicenter  to  the  ISC  cratering  epicenter.  The 
relative  locations  are  probably  only  good  to  within  a  few  km  or  so. 

Five  methods  were  used  to  estimate  the  log(Pmax/Pa)  average  for  each  event  and 
combinations  of  events.  In  the  First  method,  the  maximum  likelihood  (ML)  network 
magnitude  was  calculated  for  each  event  using  log(A)  and  log(A/T)  for  the  "a”,  "b", 
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and  "max"  phase  of  the  P  wave  at  each  station.  The  amplitude,  A,  was  corrected  for 
the  gain  of  the  WWSSN  station  at  the  estimated  period,  T.  The  network  magnitudes 
suffer  from  the  requirement  for  station  corrections.  In  the  second  method,  the 
log(Pmax/Pa)  value  was  averaged  over  the  network  for  each  event.  This  method  suffers 
from  non-detections,  clippings,  and  non-uniform  station  distributions.  In  the  third 
method,  only  common  stations  were  used  for  two  pairs  of  events.  This  method  suffers 
from  the  fact  that  due  to  station  down-time,  non-detections  and  clipping,  the  subset  of 
stations  common  to  two  events  is  always  less  than  the  number  of  stations  available  for 
each  event.  The  fourth  method,  was  a  maximum  likelihood  estimate  of  the 
log(Pmax/Pa)  over  the  network  for  each  event.  The  fifth  method,  was  a  maximum  likel¬ 
ihood  estimate  of  the  log(Pmax/Pa)  value  for  combinations  of  events.  The  fourth  and 
fifth  methods  combine  to  give  the  least  biased  estimates  of  the  value  of  log(Pmax/Pa) 
since  the  clipping  of  Pmax  and  non-dctcctions  of  Pa  may  be  accounted  for  by  the  max¬ 
imum  likelihood  procedure.  Furthermore,  the  ML  log(Pmax/Pa)  values  are  not  affected 
by  station  corrections  since  the  station  correction  is  canceled  by  the  ratio  of  Pmax/Pa. 

Network  magnitudes  were  determined  to  estimate  the  network  magnitude  differen¬ 
tials  between  mb(a),  mb(b),  and  mb(max).  The  method  of  maximum  likelihood  magni¬ 
tude  determination  (Ringdal  1976)  was  used  to  account  for  both  clipping  and  non¬ 
detections  across  the  network.  For  network  magnitude  determination,  the  Veith  and 
Clawson  (1972)  tables  were  applied  with  individual  station  corrections  from  Blandford  et 
al.  (1984).  These  station  corrections  are  listed  in  Table  Ila.  Additional  tests  were  made 
using  residuals  derived  from  a  recent  study  by  Ringdal  using  ISC  reported  magnitudes 
tabulated  in  Table  lib.  A  plot  of  the  Blandford  et  al.  (1984)  versus  Ringdal  residuals  is 
shown  in  Figure  1.  The  correlation  coefficient  between  the  two  sets  of  residuals  is  0.61. 
Because  the  two  sets  of  "station  corrections"  were  derived  in  different  ways  and  from 
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different  data  seta  we  chose  to  primarily  use  the  Blandford  et  al.  (1984)  corrections 
because  they  were  derived  from  an  explosion  data  set  with  an  effort  to  avoid  bias 
between  test  sites.  Furthermore  the  Blandford  et  al.  (1984)  station  residuals  constitute  a 
more  complete  set  of  WWSSN  stations.  Figures  2  and  3  show  the  Ringdal  (1984)  sta¬ 
tion  corrections  versus  North  (1977)  and  the  North  (1977)  versus  Blandford  et  al.  (1984) 
with  correlations  of  0.87  and  0.46  respectively.. 

Note  that  the  Ringdal  station  residuals  have  standard  errors  of  a  single  observation 
near  0.36  magnitude  units.  This  is  a  good  estimate  of  the  "typical"  width  of  the  distri¬ 
bution  of  magnitudes  from  a  network  of  stations.  Therefore,  in  order  to  see  expected 
biases  between  a  cratering  and  non-cratering  event  near  0.1  magnitude  units  with  omtm 
=  0.05  requires  on  the  order  of  40  or  more  stations  in  a  network. 

For  each  event  and  at  each  station,  the  maximum  P-wave  amplitude  in  the  first  5 
seconds  was  measured,  along  with  the  "b"  phase,  and  the  "a"  phase.  In  the  absence  of  a 
measurable  phase,  the  clipping,  or  noise  level  was  measured.  This  procedure  applied  to 
"max",  "b",  and  "a"  phase  alike.  In  general,  teleseismic  distances  between  20  and  95 
degrees  were  used,  however,  experiments  using  only  stations  30  and  90  degrees  from  fhe 
epicenter,  rarely  changed  the  network  magnitude  by  more  than  0.05  magnitude  units. 
These  distance  ranges  are  common  for  teleseismic  mb-yield  determinations.  Further¬ 
more,  log(A/T)  for  the  Pmax  phase  is  used  almost  universally  to  determine  mb.  Conse¬ 
quently  the  most  important  bias,  3,  that  can  be  determined  is  for  Pmax  using  log(A/T). 
Calculations  were  performed  for  log(A)  magnitude  as  well. 

The  network  maximum-likelihood  magnitude  results  are  shown  in  Table  Ilia,  Illb, 
and  IIIc.  The  magnitudes  were  calculated  using  both  the  log(A)  and  log(A/T)  from  the 
station  amplitude  and  period  readings. 


Teledyne  Geotech 


8 


MAY  1985 


TASK  1 


EXPLOSION  CRATERING  YIELD  DIAS 


Table  Ilia.  Maximum- Likelihood  Magnitudes 
20  <  A  <  95 

Station  Effects  From  Table  Ila. 


Event 

mb(a) 

mb(a) 

mb(b) 

mb(b) 

mb(max) 

mb(max) 

A 

A/T 

A 

A/T 

A 

A/T 

15jan65 

19jun68 

llsep69 

30nov69 

14dec73 

27apr75 

04jul76 

07dec76 

lljun78 

15sep78 


Table  Illb.  Maximum- Likelihood  Magnitudes 

30  <  A  <  90 

Station  Effects  From  Table  Ila. 

Event 

mb(a) 

mb(a) 

mb(b) 

mb(b) 

mb(max) 

mb(max) 

A 

A/T 

A 

A/T 

A 

A/T 

15jan65 

5.43 

5.50 

5.66 

5.72 

5.80 

5.87 

I9jun68 

4.56 

4.64 

5.03 

5.08 

5.23 

5.28 

llscp69 

3.71 

3.84 

4.15 

4.21 

4.49 

4.52 

30nov69 

5.36 

5.41 

5.72 

5.77 

5.90 

5.94 

14dec73 

5.16 

5.18 

5.46 

5.49 

5.70 

5.72 

27apr75 

4.76 

4.86 

5.13 

5.23 

5.27 

5.37 

04jul76 

5.15 

5.19 

5.54 

5.58 

5.86 

5.90 

07dec76 

4.85 

4.91 

5.27 

5.33 

5.44 

5.49 

lljun78 

5.13 

5.18 

5.45 

5.50 

5.73 

5.78 

15sep78 

5.35 

5.42 

5.62 

5.69 

5.78 

5.85 
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Table  111c.  Maximum-Likelihood  Magnitudes 
30  <  A  <  90 

Station  Effects  From  Table  II b. 


Event 

mb(a) 

mb(a) 

-O 

S 

E 

mb(b) 

mb(max) 

mb(max) 

A 

A/T 

A 

A/T 

A 

A/T 

15jan65 

19jun68 

llsep69 

30nov69 

14dec73 

27apr75 

04jul76 

07dec76 

lljun78 

15sep78 


The  network  magnitude  differences  for  mb(Pmax)-mb(Pa),  mb(Pb)-mb(Pa),  and 
mb(Pmax)-mb(Pb)  for  both  log(A)  and  log(A/T)  are  shown  in  Tables  IVa,  IVb,  and 
IVc.  These  tables  simply  represent  the  differences  that  can  be  calculated  from  Tables 
Ilia,  IHb,  and  IIIc,  respectively. 


Table  IVa.  Maximum-Likelihood  Differences  From  Table  Ilia. 
20  <  A  <  95 


Event 

Pmax-Pa 

Pmax-Pa 

Pb-Pa 

Pb-Pa 

Pmax-Pb 

Pmax-Pb 

(A) 

(A/T) 

(A) 

(A/T) 

(A) 

(A/T) 

19jun68 

0.65 

0.66 

0.44 

0.43 

0.21 

0.23 

IIsep69 

0.76 

0.73 

0.39 

0.38 

0.37 

0.35 

30nov69 

0.57 

0.56 

0.37 

0.36 

0.20 

0.20 

14dec73 

0.53 

0.54 

0.30 

0.31 

0.23 

0.23 

27apr75 

0.55 

0.55 

0.33 

0.33 

0.22 

0.22 

04jul76 

0.70 

0.70 

0.38 

0.38 

0.32 

0.32 

07dec76 

0.68 

0.68 

0.50 

0.50 

0.18 

0.18 

Ujun78 

0.61 

0.60 

0.32 

0.30 

0.29 

0.30 

I5sep78 

0.49 

0.48 

0.30 

0.30 

0.19 

0.18 

AVE 

0.62(.08) 

0.61(.08) 

0.37{.06) 

0.37(.06) 

0.25(.06) 

0.25(.06) 

15jan65 

0.38 

0.38 

0.23 

0.23 

0.15 

0.15 
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Table  IVb.  Maximum-Likelihood  Differences  From  Table  Illb. 

30  <  A  <  90 _ 

Event  Pmax-Pa  Pmax-Pa  Pb-Pa  Pb-Pa  Pmax-Pb 


Pmax-Pb 


AVE 

0.58(.12) 

0.56(.ll) 

0.37(.06) 

0  36(.05) 

0.23(.07) 

Emumi 

15jan65 

0.37 

0.37 

0.23 

0.22 

0.14 

0.15 

Table  IVc  Maximum  Likelihood  Magnitude  Differences  From  Table  111c. 
30  <  A  <  90 

Event 

Pmax-Pa 

Pmax-Pa 

Pb-Pa 

Pb-Pa 

Pmax-Pb 

Pmax-Pb 

(A) 

(A/T) 

(A) 

(A/T) 

(A) 

(A/T) 

19jun68 

0.66 

0.68 

0.47 

0.47 

0.19 

0.21 

llsep69 

0.77 

0.80 

0.40 

0.46 

0.37 

0.34 

30nov69 

0.56 

0.55 

0.35 

0.34 

0.21 

0.21 

14dec73 

0.53 

0.53 

0.30 

0.30 

0.23 

0.23 

27apr75 

0.53 

0.51 

0.37 

0.37 

0.16 

0.14 

04jul76 

0.70 

0.71 

0.38 

0.38 

0.32 

0.33 

07dec76 

0.57 

0.56 

0.42 

0.42 

0.15 

0.14 

lljun78 

0.61 

0.61 

0.32 

0.32 

0.29 

0.29 

15sep78 

0.43 

0.43 

0.27 

0.27 

0.16 

0.16 

AVE 

0.60(.10) 

0.60(.U) 

0.36(.06) 

0.37(.07) 

0.23(.07) 

0.23(.07) 

15jan65 

0.38 

0.37 

0.23 

0.22 

0.15 

0.15 

The  standard  deviation  ol  the  event  magnitude  difference  population  is  shown  in 
"()"  for  the  population  of  non-cratering  events.  Although  the  cratering  event  difference, 
mb(Pmax)-mb(Pa),  lies  outside  the  one  sigma  range  for  all  the  columns  in  Table  IV, 
there  are  several  events  for  which  the  relative  sizes  of  "Pa"  and  "Pmax"  are  closer  to  the 
cratering  shot  than  the  non-cratering  average.  From  the  point  of  view  of  similarity  of 
size  (mb)  and  proximity  ( 6  from  Table  I.),  if  we  take  those  events  for  which  mb  >  5.6 
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and  6  <  10  km  then  the  average  estimated  bias,  E,  between  the  the  cratering  and  non¬ 
cratering  shot  ranges  between  0.10  and  0.32  based  on  mb(Pmax)-mb(Pa)  and  0.07  and 
.27  based  on  the  mb(Pb)-mb(Pa)  values.  These  numbers  vary  depending  upon  the 
choice  of  log(A)  or  log(A/T)  or  the  comparison  of  the  "Pa"  phase  to  the  "Pb"  phase  or 
to  the  "Pmax"  phase.  Unquestionably,  the  "Pa”  phase  is,  on  the  average,  larger  with 
respect  to  the  "Pb"  or  "Pmax"  phases  for  the  cratering  event  than  for  the  general  popu¬ 
lation  of  E.  Kazakh  contained  explosions.  The  cratering  explosion  has,  on  average,  only 
a  slightly  smaller  ratio  of  of  Pmax/Pb. 

In  order  to  examine  the  relative  sizes  of  the  "Pa”  vs  "Pmax"  on  a  station-by¬ 
station  basis,  and  to  avoid  the  unnecessary  use  of  station  corrections,  the  average 
station-by-station  average  of  log(Pmax/Pa)  was  determined  for  each  event  over  the  sta¬ 
tions  detecting  both  "Pmax"  and  "Pa".  The  results  are  tabulated  in  the  following 
Table. 


Table  V.  Station  Average  n 

>b(Pmax)-mb(Pa)  With  Sigma 

Event 

log(Pmax/Pa) 

log(Pmax/Pa) 

log(A/A) 

log[(A/T)/(A)l 

Jun68 

0.43(.04) 

0.4l(.04) 

Sep69 

0.34(.05) 

NA 

Nov69 

0.52(.03) 

0.5l(.03) 

Dec73 

0.39(.04) 

0.39(.05) 

Apr75 

0.51(.05) 

Jul76 

0.58(.04) 

'-rWzL*..  IS 

Dec76 

0.56(.07) 

0.57(.07) 

Jun78 

0.50(.04) 

0.51(.04) 

Sep78 

0.42(.03) 

NA 

AVE 

0.46(.02) 

0.46(.02) 

Jan65 

0.38(.03) 

0.39(.03) 

The  standard  deviation  of  the  mean  is  shown  in  "()".  Again,  we  observe  that  the 
)og(Pmax/Pa)  is  smaller  for  the  cratering  shot.  The  difference  between  the  cratering  shot 
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and  the  non-cratering  explosion  is  £  -=■  0.08(0.04).  However  if  we  limit  the  comparison  to 
shots  with  mb  >  5.6  and  6  <  10  km  then  £  J?  0.12  magnitude  units  for  all  but  the  Scp78 
event.  The  difference  in  the  log(Pmax/Pa)  between  the  Sep78  and  Jan65  events  is  only 
0.04(0.03).  The  Sep78  event  clearly  stands  out  among  the  larger  non-cratering  events. 

Since  the  networks  may  change  in  subtle  ways  from  one  event  to  another,  we  per¬ 
formed  the  same  calculations,  network  and  station-by-station,  using  only  common  sta¬ 
tions.  Only  stations  that  were  common  to  both  the  Jan65  and  the  Dcc73  events  or  com¬ 
mon  to  both  the  Jan65  and  Sep78  events. 


Table  VI.  Network  Max-Like  Magnitude  Differences 

30  <  A  <  95 

Common  Stations  only. 

Event 

mb(max)-mb(a) 

mb(max)-mb(a) 

#  stn 

(A)fc(A) 

(A/T)*(A) 

Jan65 

0.40 

HI 

Dec73 

0.46 

mm 

Jan65 

0.33 

0.38 

22 

Sep78 

0.42 

0.51 

22 

Using  common  stations  only  for  the  Jan65-Dec73  and  Jan65-Sep78  comparisons 
increases  the  differences  between  the  log(Pmax/Pa)  averaged  over  the  entire  network. 
The  differences  are  then  0.16  and  0.09  for  log(A)  respectively.  The  Dec73  event  is  17  km 
from  the  Jan65  event  and  may  represent  a  different  testing  environment  while  the  ISC 
location  of  the  Sep78  event  is  within  10  km  of  the  Jan65  event  and  it  remains  uncertain 
why  the  event  has  a  different  relative  size  of  "Pa"  versus  "Pmax"  than  the  general  popu¬ 
lation  of  contained  explosions. 

Alternatively,  one  can  calculate  the  station-by-station  differences  for  a  common  set 
of  stations. 
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Table  VII.  St&tion-by-Station  Magnitude  Differences 

30  <  A  <  95 

Common  Stations  Only 


Event 

log(Pmax/Pa) 

log(Pmax/Pa) 

#  stn 

log(A/A) 

log[(A/T)/(A)] 

Jan65 

0.37(.04) 

0.42(.04) 

28 

Dec73 

0.46(.04) 

0.50(.04) 

28 

Jan  65 

0.18(.03) 

0.38(.05) 

25 

Sep78 

0.42(.04) 

0.53(.06) 

25 

The  standard  deviations  of  the  mean  arc  shown  in  Using  the  data  of  Table  VII, 

above,  the  differences  between  the  Jan65  and  Dec73  and  Jan65  and  Sep78  events  for 
log(Pmax/Pa)  are  0.09(.04)  and  0.22(.04)  respectively.  The  differences  expressed  in 
Tables  VI  and  VII  above  are  greater  than  expected  given  the  quoted  standard  deviations 
of  Table  VII. 

In  order  to  eliminate  the  need  for  station  corrections,  and  to  account  for  the  biasing 
effects  of  non-detection  of  Pa,  and  clipping  of  Pmax,  the  maximum-likelihood  estimation 
technique  was  used  to  estimate  the  log(Pmax/Pa)  over  a  network  of  stations.  Several 
cases  could  be  added  to  the  case  where  measurement  of  Pmax  and  Pa  could  be  made. 
This  method  has  the  advantage  that  there  is  a  greater  probability  that  the  "Pa"  phase 
may  not  be  observed  while  the  "Pmax"  phase  is  observed.  Consequently,  if  the  correla¬ 
tion  of  "Pa"  and  "Pmax"  amplitudes  is  not  100%  then  there  will  be  preferential  sam¬ 
pling  of  stations  where  the  ratio  (Pmax/Pa)  is  smaller  than  the  norm.  These  cases  are 
summarized  in  the  following  table. 
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Table  VIII. 

Pmax  on  scale 

Pmax  clipped 

j  Pmax  <  noise 

Pa  on  scale 

signal 

clipped 

not  used 

Pa  clipped 

not  used 

not  used 

not  used 

Pa  <  noise 

clipped 

clipped 

not  used 

Figures  4,  5  and  6  illustrate  the  use  of  the  maximum  likelihood  approach  to  the 
estimation  of  the  statistic  log(Pmax/Pa).  The  arrows  on  the  histograms  represent  the 
number  of  clipped  values  in  the  population  of  log(Pmax/Pa)  measurements.  The  follow¬ 
ing  results  were  found  for  the  ten  individual  events. 


Table  IX.  Maximum- 

Likelihood  log(Pmax/Pa)  Ratios 

Event 

log(Pmax/Pa) 

19jun68 

0.68(0.33) 

llsep69 

0.77(0.34) 

30nov69 

0.54(0.20) 

14dec73 

0.51(0.25) 

27apr75 

0.60(0.25) 

04jul76 

0.72(0.25) 

07dec76 

0.89(0.44) 

lljun78 

0.53(0.21) 

15sep78 

0.47(0.28) 

AVE 

0.63(0.13) 

15jan65 

0.41(0.24) 

The  "()”  standard  error  given  for  each  event  is  the  maximum  likelihood  estimated 
population  sigma  for  the  log(Pmax/Pa)  over  the  network.  The  standard  error  given  for 
the  "AVE”  in  the  table  is  the  standard  error  for  the  non-cratering  explosion  population. 
These  results  are  in  general  agreement  with  the  results  from  Tables  V,  VI,  and  VII.  The 
average  bias  is  estimated  as,  5  =  0.63  -  0.41  =  0.22,  using  all  nine  contained  explosions. 
Figures  7,  8,  and  9  compare  the  composite  maximum  likelihood  fits  to  the  combined 
populations  of  7  non-cratering  events  5  non-cratering  events  and  4  non-cratering  events 
respectively.  The  7  largest  non-cratering  events  have  a  combined  population  of 
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log(Pmax/Pa)  with  an  estimated  mean  of  0.58  (Figure  7)  which  implies  E  —  0.17.  The  5 
non-cratering  events  (30  Nov  69,  14  Dee  73,  04  Jul  76,  07  Dee  76,  and  15  Sep  78)  have  a 
combined  population  of  log(Pmax/Pa)  with  an  estimated  mean  of  0.60  (Figure  8)  which 
implies  E  =  0.19.  The  results  for  four  events  are  shown  in  Figure  9  with  an  estimated 
mean  of  0.68  and  an  implied  E  =  0.27.  The  four  non-cratering  events  chosen  for  Figure 
9  are  the  same  events  used  by  Dcr  et  al.  (1985)  in  for  source  deconvolution  (Task  2,  this 
contract)  using  the  Shumway-Der  multichannel  deconvolution  technique.  These  four 
events  are  lljun78,  7jul76,  4jul76,  and  27apr75.  This  analysis  demonstrates  the  variabil¬ 
ity  of  the  statistic  log(Pmax/Pa)  and  implies  that  E  lies  in  the  range  between  0.17  and 
0.27. 

Of  interest  is  that  the  typical  population  width,  a,  —  0.29  ,  for  log(Pmax/Pa)  is 
only  20%  less  than  the  population  width,  a0  =  0.36,  of  station  magnitude  residuals.  The 
average  station  residual  a  from  Table  lib  is  about  0.36  magnitude  units.  The  individual 
station  values  of  log(Pmax/Pa)  contain  a  great  deal  of  remaining  scatter;  although  the 
amplitude  of  Pa  is  correlated  to  the  amplitude  of  Pmax,  the  correlation  is  imperfect. 

CONCLUSIONS 

In  conclusion,  it  appears  that  log(Pmax/Pa)  is  variable  from  shot  to  shot  but  that 
it  generally  is  between  0.5  and  0.7  magnitude  units.  The  P  waveforms  from  the  crater¬ 
ing  event  of  Jan  15,  1965,  on  average,  have  values  of  log(Pmax/Pa)  generally  near  to  0.4. 
The  estimated  bias  to  correct  the  mb(Pmax)  for  the  cratering  event  to  a  contained  shot 
with  a  yield  of  125  Kt  is  then  E  =0.10  to  .30  magnitude  units.  The  preferred  values 
depending  upon  the  method  used  to  estimate  the  difference  E  =  log  (Pmaz  /  Pa)„, 

•  log  (Pmax  /  Pa  )t,««(r .  The  maximum  likelihood  estimates  for  the  individual  network 
values  of  log(Pmax/Pa)  from  Table  IX  are  the  most  direct  and  least  biased  estimates  of 
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the  log(Pmax/Pa)  values  for  each  event.  These  estimates  are  the  least  biased  by  non- 
detection  of  Pa  and  clipping  of  Pmax.  The  estimates  of  Table  IX  are  non  influenced  by 
the  choice  of  station  corrections,  and  they  are  not  influenced  by  the  choice  of  Veilh  and 
Clawson  distance  corrections. 

The  maximum  likelihood  estimates  of  log(Pmax/Pa)  range  from  0.58  to  0.68 
depending  upon  the  set  of  contained  explosions  that  are  chosen.  The  best  estimate  of 
the  network  mean  of  the  log(Pmax/Pa)  for  the  cratering  explosion  is  0.41.  From  these 
values,  the  estimated  mb(Pmax)  bias  for  the  cratering  explosion,  E,  ranges  from  0.17  to 
0.27.  From  Table  Ilia  the  mb(Pmax)  for  the  cratering  event  is  5.88,  and  the  implied 
unbiased  mb(Pmax)  estimate  of  the  15  Jan  65  cratering  event  is  then  6.05  to  6.15.  If  we 
assume  that  the  cratering  explosion  had  a  yield  of  125  Kt,  then  a  log(yicld)-mb  slope  of 
1.0  would  give  an  estimate  of  6.13  to  6.23  for  the  150  Kt  scale-depth  shot.  Given  the 
uncertainties  in  the  magnitude  determinations  and  the  uncertainties  in  the  scaling  laws 
for  P-fpP  interference  the  conclusions  are  consistent  with  the  estimated  mb(150  Kt) 
value  of  6.17  given  by  Blandford  et  al.  (1984).  It  should  be  emphasized  that  such  a 
prediction  is  based  on  the  WWSSN  network  using  distances  between  20  and  90  degrees, 
the  yield  value  of  125  Kt  for  the  15  Jan  65  cratering  explosion,  and  the  assumption  that 
the  Pa  phase  amplitude  is  the  same  for  the  cratering  event  and  a  contained  explosion  of 
the  same  yield. 

However,  events  can  be  found  within  the  East  Kazakh  test  site  with  anomalous 
log(Pmax/Pa)  values  on  the  average  over  a  network  of  teleseismic  stations.  The  event  of 
15  Sep  78  had  a  network  average  value  of  !og(Pmax/Pa)  that  was  only  0.06  (from  Table 
IX)  larger  than  the  cratering  event.  This  investigation  raises  the  question  as  to  what 
mechanism  could  be  responsible  for  the  reduction  in  the  Pmax  of  this  event  (or  alterna- 
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tivcly  larger  Pa).  Perhaps  other  information  is  available  that  indicates  whether  the  15 
Sep  78  event  was  unusual.  A  cursory  investigation  of  the  WWSSN  waveforms  the  the 
East  Kazakh  events  studied  did  not  reveal  a  marked  difference  between  the  waveforms  of 
the  15  Sep  78  and  other  events.  Furthermore,  it  is  not  always  clear  whether  the  crater¬ 
ing  event  is  different  from  the  remaining  population  of  events  until  statistical  averaging 
is  performed  on  the  log(Pmax/Pa)  statistics;  such  is  the  stochastic  nature  of  short  period 
waveforms. 

North  and  Fitch  (1981)  analyzed  the  surface  waves  of  the  15  Sep  78  and  11  June  78 
events  as  well  as  several  others.  They  found  the  15  Sep  78  event  to  be  of  "TYPE  1";  the 
Rayleigh  waves  were  not  reversed  at  any  observed  azimuth.  According  to  the  North  and 
Fitch  (1981)  epicenter  relocations,  the  11  June  78  and  the  15  Sep  78  events  are  closer  to 
each  other  than  each  is  to  the  cratering  event  epicenter  (see  Figure  10  taken  from  North 
and  Fitch,  1981).  The  11  June  78  event  appears  to  have  a  normal  log(Pmax/Pa)  distri¬ 
bution  over  the  network. 

Herrin  and  Goforth  (1982)  analyzed  the  Rayleigh  waves  of  the  15  sep  78  event  and 
although  they  did  not  observe  reversed  Rayleigh  waves  for  the  15  Sep  78  explosion,  they 
observed  larger  arrivals  proceeding  the  main  fundamental  at  KAAO  and  an  apparent 
delay  of  the  Rayleigh  wave  by  3  to  4  seconds  at  SHIO.  How  this  applies  to  the 
log(Pmax/Pa)  statistics  is  uncertain,  for  it  is  not  clear  how  the  Pmax  amplitude  would 
be  relatively  smaller  for  an  event  with  tectonic  complications. 
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aae 

-0.34 

aam 

0.24 

ade 

0.22 

aft 

-0.10 

aku 

-0.10 

alq 

-0.25 

anp 

-0.56 

ant 

0.05 

aqu 

-0.24 

are 

0.11 

atl 

0.14 

atu 

-0.12 

bag 

-0.09 

bdf 

-0.01 

bee 

-0.12 

bhp 

-0.05 

bks 

0.14 

bla 

0.18 

bog 

-0.02 

boz 

0.07 

bul 

0.08 

car 

0.08 

chg 

-0.23 

cmc 

-0.31 

col 

0.21 

cop 

0.02 

cor 

0.27 

eta 

0.01 

dag 

0.00 

dal 

0.30 

dav 

-0.14 

dug 

0.05 

eil 

0.01 

eak 

-0.03 

Ho 

0.39 

fvm 

0.19 

gdh 

-0.07 

geo 

-0.17 

gie 

-0.24 

gol 

-0.17 

grm 

-0.24 

gsc 

-0.02 

gua 

0.13 

hkc 

-0.23 

hlw 

-0.11 

hn- 

0.21 

hnr 

0.22 

ist 

0.08 

jet 

0.10 

jer 

-0.07 

kbl 

-0.02 

kba 

-0.01 

kev 

-0.07 

kip 

0.20 

kod 

0.12 

kon 

0.13 

krk 

-0.27 

ktg 

-0.18 

lah 

0.25 

lem 

-0.29 

Ion 

-0.09 

lor 

0.13 

lpa 

0.05 

Ipb 

-0.13 

Ips 

-0.14 

lub 

0.39 

mal 

-0.05 

man 

0.59 

mat 

-0.01 

mds 

0.27 
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FIGURE  CAPTIONS 

FIGURE  1.  Blandford  at  al.  (1984)  stations  corrections  versus  Ringdal  (1984)  station 
corrections.  Correlation  coefficient  of  0.65. 

FIGURE  2.  North  (1977)  station  corrections  versus  Ringdal  (1984)  station  corrections. 
Correlation  coefficient  of  0.87. 

FIGURE  3.  North  (1977)  station  corrections  versus  Blandford  et  al.  (1984)  station 
corrections.  Correlation  coefficient  of  0.46. 

FIGURE  4.  Histogram  of  the  log(Pmax/Pa)  measurements  for  the  30  Nov  69  event. 
Clipped  values  of  the  log(Pmax/Pa)  are  indicated  by  right  pointing  arrows.  There  are  20 
signal  measurements,  and  10  clipped  measurements.  The  maximum  likelihood  estimate 
of  the  Gaussian  distribution  is  plotted  with  a  mean  of  0.54  and  a  a  —  0.20. 

FIGURE  5.  Histogram  of  the  log(Pmax/Pa)  measurements  for  the  14  Dec  73  event. 
Clipped  values  of  the  log(Pmax/Pa)  are  indicated  by  right  pointing  arrows.  There  are  27 
signal  measurements,  and  23  clipped  measurements.  The  maximum  likelihood  estimate 
of  the  Gaussian  distribution  is  plotted  with  a  mean  of  0.51  and  a  o  =  0.25. 

FIGURE  6.  Histogram  of  the  log(Pmax/Pa)  measurements  for  the  15  Sep  78  event. 
Clipped  values  of  the  log(Pmax/Pa)  are  indicated  by  right  pointing  arrows.  There  are  21 
signal  measurements,  and  10  clipped  measurements.  The  maximum  likelihood  estimate 
of  the  Gaussian  distribution  is  plotted  with  a  mean  of  0.47  and  a  a  =  0.28. 

FIGURE  7.  Composite  histogram  for  the  7  largest  non-cratering  events.  106  signals  and 
69  clipped  measurements.  Maximum  likelihood  estimate  of  the  Gaussian  distribution 
has  mean  0.58  and  a  -*  0.29.  Distribution  for  the  15  Jan  65  cratering  event  is  shown 
above  with  mean  0.41  and  a  —  0.24. 
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FIGURE  8.  Composite  histogram  for  the  5  largest  non-cratering  events.  72  signals  and 
55  clipped  measurements.  Maximum  likelihood  estimate  of  the  Gaussian  distribution 
has  mean  0.60  and  a  =  0.29.  Distribution  for  the  15  Jan  65  cratering  event  is  shown 
above  with  mean  0.41  and  a  =  0.24. 

FIGURE  9.  Composite  histogram  for  4  non-cratering  events.  38  signals  and  26  clipped 
measurements.  Maximum  likelihood  estimate  of  the  Gaussian  distribution  has  mean 
0.68  and  a  =  0.33.  Distribution  for  the  15  Jan  65  cratering  event  is  shown  above  with 
mean  0.41  and  a  =  0.24. 

FIGURE  10.  Event  relocations  of  selected  E.  Kazakh  explosions  taken  directly  from 
North  and  Fitch  (1981).  Crosses  indicate  the  relocation  from  ISC  locations  indicated  by 
small  circles.  Events  6/11/78,  9/15/78,  and  1/15/65  are  labeled  A,  B,  and  C  respec¬ 
tively. 
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